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We show that the intermediate scattering function S„ (q, t) for neutron scattering (ns) can be 
expanded naturely with respect to a set of molecular correlation functions that give a complete 
description of the translational and orientational two-point correlations in the liquid. The general 
properties of this expansion are discussed with special focus on the g-dependence and hints for a 
(partial) determination of the molecular correlation functions from neutron scattering results are 
given. The resulting representation of the static structure factor S n (q) is studied in detail for 
a model system using data from a molecular dynamics simulation of a supercooled liquid of rigid 
diatomic molecules. The comparison between the exact result for S n {q) and different approximations 
that result from a truncation of the series representation demonstrates its good convergence for the 
given model system. On the other hand it shows explicitly that the coupling between translational 
(TDOF) and orientational degrees of freedom (ODOF) of each molecule and rotational motion of 
different molecules can not be neglected in the supercooled regime. Further we report the existence 
of a prepeak in the ns-static structure factor of the examined fragile glassformer, demonstrating that 
prepeaks can occur even in the most simple molecular liquids. Besides examining the dependence 
of the prepeak on the scattering length and the temperature we use the expansion of S n (q) into 
molecular correlation functions to point out intermediate range orientational order as its principle 
origin. 

PACS numbers: 61.12.-q, 61.20.-p, 61.25.Em 
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I. INTRODUCTION 

Neutron scattering is one of the most important tools 
to determine the structure and dynamics of condensed 
matter. One of its major advantages is that the neutron 
as an electrically neutral probe is not influenced by the 
electron cloud of the target atoms but interacts only with 
the nucleus. Consequently neutron scattering can be di- 
rectly interpreted in terms of the atomic structure and 
motion. However, if one wants to examine a molecular 
system not only with respect to the constituent atoms 
but in terms of the molecular units some care has to be 
taken. Since the superposition of scattering from atomic 
sites is in general not equivalent to scattering from the 
molecular center of mass, neutron scattering from molec- 
ular liquids is sensitive to orientational as well as trans- 
lational correlations in the system. 

This was realized already many years ago p] pj and 
led to different formulations of neutron scattering by 
molecules. For reasons of simplicity and also because 
the center of mass correlations were of principle interest 
those approaches tried to " correct" for the effect of ori- 
entational degrees of freedom. In the work of Sears ||] 
this is done using the so called weak hindering approxi- 
mation, i.e. the assumptions of statistical independence 
of i) rotational and translational motion of any molecule 
and ii) rotational motion of any two different molecules. 

In recent years the mode coupling theory of the glass 
transition (see e.g. P-pJ) stimulated a lot of experiments 
on supercooled liquids. Among other techniques also 
neutron scattering was used (see e.g. B and for more re- 



cent work pftp) ) to study the static and dynamic prop- 
erties of glassforming substances most of which are of 
molecular nature, like e.g. glycerol, salol or orthoter- 
phenyl. In the situation of the dense packed molecular 
liquid where the motion becomes more and more coop- 
erative the assumptions of Sears are not reasonable any- 
more as was shown by Chen et al. |]12| for incoherent 
neutron scattering in simulated supercooled water. Be- 
low, in the study of a liquid of diatomic rigid molecules 
we will come to the same conclusion for coherent scat- 
tering. This raises the question how neutron scattering 
data can be analyzed in terms of the molecules in this 
case. 

With increasing computer power the simulations of su- 
percooled liquids trend to deal with more realistic model 
systems @ like SPC/E water or orthoterphenyl [JTJ). 
Thus it has become feasible to examine in detail also the 
orientational degrees of freedom in the strongly super- 
cooled regime (see e.g. 0,0). In experiments it is up 
to now only possible to measure orientational correlators 
for q = 0, i.e. one does not get information about the spa- 
tial arrangement of orientationally correlated molecules. 
On the other hand this information would be valuable 
to understand the interaction between translational and 
orientational degrees of freedom. Therefore it seems sen- 
sible to examine in which way orientational correlations 
enter into the neutron scattering function. 

A feature of special interest in supercooled liquids and 
glasses is the appearance of a prepeak in the static ns- 
structure factor at a g-value that corresponds to dis- 
tances larger than the average nearest neighbour dis- 



tance. As a sign of intermediate range order prepeaks 
have been studied in a variety of systems (see e.g. 19 2lJ 
and references therein). Mostly they have been at- 
tributed to the network structure of strong glasses or 
charge effects in ionic glassformcrs. The behaviour of pre- 
peaks in molecular systems seems to be quite different as 
was pointed out recently by Morineau et al. p2| in their 
study of the temperature- and pressure-dependence for 
m-toluidine and m-fluoroaniline which will be discussed 
in section 
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We will present results for the ns-static structure fac- 
tor of a liquid of diatomic molecules, i.e. one of the most 
simple systems in which effects of the molecular nature 
can be studied. We will discuss the dependence of the 
prepeak on temperature and scattering length and make 
comparison to findings for other systems. Using the rep- 
resentation of S n (q) in terms of the molecular correlation 
functions we will also give strong evidence that the ex- 
istence of the peak is closely connected to oricntational 
degrees of freedom. 

The paper is organized as follows. In Section || we 
will summarize shortly the approach of Sears to the in- 
termediate scattering function S n (g, t) for molecular sys- 
tems. We will show that dropping the weak hindering 
approximation leads naturely to a description in terms 
of molecular correlation functions. Some aspects of the 
representation for arbitrary molecules are discussed and 
hints for an application to experiments are given. Sec- 
tion III presents the results from the molecular dynamics 
(MD) simulation of a system of diatomic molecules. It 
is divided into three subsection concerning respectively 
the molecular correlation functions, the ns-static struc- 
ture factor and a detailed discussion of the prepeak. The 
final section gives some conclusions. 



II. NEUTRON SCATTERING FUNCTION 

We consider a set of N rigid molecules of identical ge- 
ometry each consisting of s atoms. The assumption of 
rigidity is one that is commonly used as well in the the- 
oretical analysis as in computer simulation of the liq- 
uid state. It is justified on the ground that at the tem- 
perature of the liquid only the lowest vibrational states 
are populated. The starting point of an analysis of neu- 
tron scattering from molecules is the following site-site- 
representation of the intermediate scattering function ( 
see e.g. |@) 

N s 

S„(q,t) = — J2 J2 ( a ioh a ioh + ( a inc) 2 5n'$vv 



i?-(4v(t)--R,-'„')\ 



(1) 



where a?" oh (aj nc ) is the coherent (incoherent) scattering 
length and Rj V {t) the position of atom v in molecule 



j at time t. The brackets (...) denote the canonical 
averaging over initial conditions. Unlike the geometry 
which is supposed to be identical for all molecules the 
scattering length are allowed to differ from molecule to 
molecule. This assumption is quite realistic since the 
chemical structure is independent of the isotopic compo- 
sition whereas the scattering lengths are. 

Appealing as this site-site-description is, it is not 
very favourable from a molecular view. The rigid 
molecules are most suitably described by the center of 
mass coordinates Rj(t) and the Euler angles £lj(t) — 
(4>j(t),9j(t),Xj(t)) giving the orientation of the body 
fixed system relative to the laboratory frame. 

To make a connection between these two different 
views Sears [0 introduced center of mass and relative 
coordinates by Rj V = Rj + rj V and used the Rayleigh 
expansion of the plane waves 



Jq-r^(t) 



= 53[M« + i)]**'ii(^) 



In 
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(2) 



where (r v ,D v ,ip v ) are the polar coordinates of atom v 
with respect to the body fixed frame with origin at the 
center of mass of molecule j. Without loss of generality 
the z-axis of the laboratory frame has been chosen to 
point in direction of q. The special functions appearing 
in eq. (0) are the spherical Bessel functions ji , the spheri- 
cal harmonics Y™ and the Wigner functions T> l mn p4,p5[ . 
In the summation of eq.(0) I runs over 0, 1,2, ... and n is 
restricted by — I < n < I. 

Following this strategy and not doing any factoriza- 
tion of translational and rotational motion or of rota- 
tional motion of different molecules one is naturely led 
to consider correlation functions of the following kind 

Sta,^(fl,*)=< I '-'[(2i + l)(2l , + l)]'^X; • 

n' 

■ (e-^V-^VlMt))® 1 ^,^,)), (3) 

which equal jf(pt mn (q,t) pi> mn <(q,0)}, the correlation 
function of the tensorial one particle density mode 
Pimn{q,t) |26|. As usual it is possible to separate them 
into self (j = j 1 ) and distinct terms (j ^ j') giving 



Sln,Vn'\Q't) — SlnMn'il't) + S ln ;/„/(<?, t). 



(4) 



This set of correlation functions is the generalization to 
arbitrary molecules of the correlators introduced in the 
mode coupling theory for a single molecule in an atomic 
liquid |27| and for liquids of linear molecules [E6[ and 
has been used recently to study molecular correlations in 
supercooled water p8| and diatomic molecules Cq] . Sim- 
ilar correlation functions have also been used earlier in 
the study of molecular liquids fl23|,p5| and go back to the 



work on statistical mechanics of nonspherical molecules 
by Steele ||. 

On first sight one might be scared by the large number 
of indices and not willing to consider an infinite set of 
correlation functions. Therefore we want to stress that 
due to the tensorial nature of the quantities (||) the sym- 
metry of the molecule may reduce the number of corre- 
lators. Secondly one can expect that reasonable results 
can be achieved already if one considers only the finite 
number of correlators with 1,1' < l co for small l co . For 
a system of diatomic molecules we will see in section [II 
that l co = 4 turns out to be sufficient to describe S n (q) in 
the g-range of interest. Depending on the size and form 
of the molecules it might be necessary to consider also 
larger values for l co . This becomes clear if one looks at 
the angular dependence of the Wigner functions since ex- 
tending the range of I corresponds mainly to improving 
the angular resolution. In this argument we have con- 
sidered the factors T> l mn (Q) in eq.(g) as "weights" which 
select molecules within a certain range of orientations for 
the canonical averaging. While for instance Vq (H) = 1 
gives equal weight to all orientations, 2?Q (f2) oc Y®(9, <f>) 
puts emphasis on molecules with 9 rs or 9 w n. This 
may elucidate a bit more the meaning of the g-dependent 
orientational correlators. 

Having introduced the suitable correlation functions 
cqs.(Q)-(^) can be easily combined to give 



Sn(q,t) =^^K: w ( 9 )si:;,(9,i) 



(d)0 



^(J)C(«'*)1 

with the coefficients 

1 Q { J 



(5) 



binMw (q) = -Y) 4ttJi (qr v )j v {qr v , ) 



Yr^vJYFV'MWcohO&h- 



(7) 



Here x — l/N^j & denotes the average over molecules. 
For the special case of a single linear molecule in an 
isotropic liquid a similar result has been given by Fra- 
nosch et al. [f27f . 

The representation eq. rtq) has a number of interesting 
properties. Immediately obvious is that it can be Fourier 
transformed to give a corresponding relation for the spec- 
tra or susceptibilities with the same coefficients. Further 
we notice that only correlators with m = enter into 
S n (q, t) which is a consequence of the isotropy of the fluid. 
From the general expression one can also draw some con- 
clusions about the g-dependence of the ns-intermediate 
scattering function. Eq.(|J) clearly reflects that neutron 
scattering from molecular liquids is the superposition of 



two contributions: i) the correlations within a molecule 
which are expressed in the incoherent terms involving the 
self-part of the molecular correlation functions and ii) 
the correlation between molecules which enter the dis- 
tinct terms. The q-dependence of both contributions 
again has two sources. The coefficients &/™ c ;/„'(g) and 
bf° h ii n ,(q) which could be termed incoherent and coher- 
ent molecular form factors respectively are completely 
determined by the molecular geometry and the scatter- 
ing lengths. Thus they can be calculated exactly if the 
molecular units are known. Their q-dcpcndcnce is given 
by the spherical Bessel functions ji(qr v ) and thus con- 
nected primarily to the distances r v of the atoms from 
the center of mass of the molecule. We will discuss the q- 
dependence of the molecular form factors in more detail 
in the following section. The more important quantities 
that enter into the (/-dependence of S n (q, t) are of course 
the molecular correlation functions which give a statis- 
tical description of the interactions and the dynamic of 
the system. We will discuss their q-dependence for the 
mentioned model system in the first subsection of the 
following section. 

The coefficients bf° h lln ,(q) and b\™ c v n , (q) represent 
weighting functions that determine at which g-values a 
given molecular correlation function Si n Vn i (q, t) appre- 
ciably contributes to S n (q,t). In turn knowledge of the 
molecular form factors could be used to attribute the 
structure of the ns-intermediate scattering function to 
some molecular correlation functions. To undertake such 
an effort one would have to do a series of experiments in 
which the molecular form factors bf^\, n , (q) and b % {^ Vn , (q) 
are varied systematically. A common technique to do so 
would be the use of isotopes of different scattering length, 
i.e. to combine results obtained with mixtures of differ- 
ent isotopic composition. From the discussion of eq.(|5|) 
given above we can conclude the following limitations: 
i) only molecular correlation functions Sf n i> n /{q,t) with 
m = are accessible , ii) information about the cor- 
relator Si nl , n ,(q,t) can only be extracted from S n (q,t) 
in a g-range where the form factor bi n j,i n i (q) is different 
from zero and not too small compared to the other form 
factors. Upon changing the scattering lengths one can 
separate the site-site correlation functions contributing 
to S n (q,t). Since these are connected by a linear rela- 
tion with the functions Si nt i' n i(q) (as a special case of 
eqs.dq)-(M)) we have Hi) the number of molecular corre- 
lation functions that can be determined (for a given q) is 
restricted by the number of site-site correlators. 

Since the molecular form factors depend also on the 
atomic configuration within the molecule another idea 
would be to compare results for molecules with similar 
geometry. To examine the structure of a liquid of di- 
atomic molecules one could for example try to combine 
results for F2, CI2, Br^ and I2 taking into account the 
shift in the average nearest neighbour distance due to the 



difference in atomic size by rescaling the g- values accord- 
ingly. In such a way one could change systematically the 
intramolecular distance d. Although this would offer a 
way to bypass the limitations given by Hi) changing the 
geometry will also affect the molecular correlation func- 
tions themselves. 

In spite of this limitations we think that experimental 
efforts in the direction of determining molecular corre- 
lation functions from neutron scattering are worth con- 
sidering. Here one has to keep in mind that up to now 
information about the orientational degrees of freedom is 
restricted to the few q = O-correlation functions that can 
be measured by dielectric response or NMR. Especially to 
understand the interaction between rotational and trans- 
lational degrees of freedom the g-dependent molecular 
correlation function have to be considered. 



III. RESULTS 

We will now turn to the examination of a special sys- 
tem. On the one hand this will give us the opportunity to 
illustrate the molecular correlation functions and the rep- 
resentation of the scattering function introduced in sec. 
|J. On the other hand the system shows interesting fea- 
tures which are apparent in the ns-static structure factor 
S n (q) and give on their own a motivation for a detailed 
inspection. 

We will study a simulated supercooled liquid of 500 
diatomic rigid molecules each consisting of two atoms 
which will be labeled A and B. The atoms have equal 
masses but the head-tail-symmetry is slightly broken by 
the interactions which are given by a superposition of 
Lennard-Jones potentials between the atomic sites. The 
molecular bond length was fixed at d — 0.5 in units of 
the Lennard-Jones radius of the A-atoms which we will 
use throughout this article. Particularly, q is given in 
units of 2-7T times the inverse of that radius. For further 
details about the molecular dynamics simulation and the 
potential we refer the reader to refs. pq| . 

In the following we will concentrate on the static prop- 
erties of the system. This is done because of a number 
of reasons: i) To test the convergence of the series rep- 
resentation (||) as function of the /-cutoff, it is sensible 
to consider the worst case. Usually this should be given 
by the static case as the following argument explains. 
If one assumes that the relaxation of the ODOF takes 
place primarily through small angular variations and that 
large angular jumps can be neglected the correlators with 
larger I, i.e. better angular resolution, will decay faster 
than those with small I. Therefore one concludes that 
for t > terms with larger I are less important than in 
the static case, ii) Last but not least we are primarily 
interested in the q dependence. Thus it will be necessary 
to choose a fixed time t and t = being a natural choice. 



The readers interested as well in the dynamic of the 
studied system and especially experimentators who are 
interested in the timescales involved shall be referred to 
the publications of Kammerer et al. flfl] where these as- 
pects have been discussed for TDOF as well as ODOF. 



A. Molecular Correlation Functions 

We will now present molecular correlation functions for 
the system of diatomic rigid molecules. As already men- 
tioned in section |nj the number of independent correlators 
can be reduced for molecules possessing an intrinsic sym- 
metry. In our case we deal with linear molecules that are 
invariant under any rotation around the axis connecting 
atoms A and B. This axis will be chosen as the z-axis of 
the body fixed frame of reference in the following. 

Symmetry considerations similar to those carried out 
in ref. J28| show that the distinct part of the static cor- 
relation functions fulfills 



InA'n 



'(?) 



S$ m {q) 



n = ri = 
otherwise 



(8) 



Since the self part of the static correlation is given sim- 
ply by Sfayni (q) = Sw 5 nn i all the important structure is 
contained in the correlation functions with n = n' = 0. 
As pointed out in section O for neutron scattering only 
correlators with m = are relevant. Thus in the fol- 
lowing we will only consider the quantities Sw (q) = 
S? [i (q)- We can further restrict ourselves to I < I' since 
the correlation functions Sw (?) are real and symmetric 
with respect to I and I' as can be easily seen from their 
definition (g). 

The correlation functions up to I, I' < l co — 2 have al- 
ready been presented by Kammerer et al. pq| . We want 
to take the opportunity to point out a minor error in this 
publication: In order to get the correct data which is in 
accordance with the chosen conventions (ours and theirs) 
the off diagonal terms given in Fig. 3 of the last reference 
of |16[ have to be multiplied by -1. The correct graphs 
for the m = 0-terms are given in our Figjj] (b). 

In addition we have determined all correlators up to 
1,1' < lco = 4. Figures fi] (c)-(e) show the m = 0-terms 
which are relevant for neutron scattering. In the discus- 
sion given below we will only refer to this data although 
the properties of the m ^ 0-terms (not shown here) are 
quite similar. 

All the molecular correlation functions have been eval- 
uated at T = 0.477 (in units of the Lennard-Jones energy 
caa of the A-particles) which is the lowest temperature 
considered in the simulation. This is a temperature lo- 
cated in the supercooled regime very close to the critical 
temperature T c = 0.475 of the mode coupling theory pw. 

We will now give a discussion of the molecular corre- 
lators Sw (?) by comparing their g-dependence. Thereby 



we want to point out the similarities and differences be- 
tween the g-dependent orientational correlators (I and/or 
V different from zero) and the better known center of 
mass correlations (/ = V = 0). We also mention some 
features for which a general relation between the q- 
dependence and the values of I and I' seems to exist. 
Thereby we will substantially enlarge the discussion given 
by Kammerer et al. which focussed mainly on the in- 
dependence and the effects of the approximate head-tail 
symmetry of the molecule. 

The most prominent peak is displayed by the cen- 
ter of mass correlator Soo(q). Located at a g-value of 
Imax ~ 6.6 it represents the first order of the nearest 
neighbour peak. Such a peak can also be found at al- 
most exactly the same position for all correlation func- 
tions Soi'(q) with I = 0. Also for the other correlators 
this nearest neighbour peak exists but it is usually shifted 
to a slightly different value than q m ax- This shift can 
easily be understood. Whereas for Soo{q) an average 
over all possible orientations of the molecule is done the 
factor Doq^j^oq^j') in eq.([|) for I and/or I' unequal 
zero "restricts" the average to a certain range of orienta- 
tions. Since in the supercooled liquid the molecules are 
very closely packed the characteristic distance between 
molecules depends on the choice of orientations Qj and 
flji simply because of steric hindrance. Thus for different 
correlators a shift in the position of the peaks is expected. 
A strong support for this picture is also that the peak is 
much less shifted for Soi'(q) as for Sw(q) with I and /' 
unequal zero because an average over all possible relative 
orientations is already done if one of the angles Qj or Qji 
is not restricted. 

Besides the shift in the peak position we observe that 
the peak amplitude also depends strongly on the cor- 
relator. One reason already pointed out by Kammerer 
et al. fllql is the approximate head-tail-symmetry of the 
molecule which results in a smaller amplitude for cor- 
relators with I + I' odd. Further we can state that the 
amplitude of the main peak tends to decrease with in- 
creasing 1,1'. 

Apart from the main maximum the center of mass cor- 
relation function Soo{q) shows also structure at higher 
values of q. Clearly perceptible is a double peak at 
q w 10.5 and q « 13.2. Whereas the peak at 5 « 13.2 
can be identified as the second order of the nearest neigh- 
bour peak the origin of the peak at q sa 10.5 is not quite 
clear. It could also be the second order of a peak located 
around q sa 5.25 which is merged with the main peak. 

Rich structure at higher values of q is also found in all 
other molecular correlation functions. If one takes into 
account possible shifts in the peak position the peaks lo- 
cated at q ~ 12 — 14, q rs 16 — 19 and q w 23 — 25 can 
be attributed to the higher orders of the nearest neigh- 
bour peak. While their principle origin thus seems to 
be clear the interesting structure that appears in their 
shape, which in some cases clearly indicates a double 



peak, can not be understood from the present investiga- 
tion. 

A further aspect which immediately strikes the eye 
is that, while the oscillations at higher q are strongly 
damped for small I and I', the higher order peaks are 
comparable or even larger than the first maximum for 
the correlators with I, I' G {2, 3, 4}. 

So far we have only considered the structure of the 
molecular correlation functions for g-values larger than 
the value of the first maximum corresponding to dis- 
tances smaller than the average nearest neighbour dis- 
tance. The center of mass correlator Soo(q) for q < 6.6 re- 
sembles the structure also found for simple liquids which 
is in strong contrast to the behaviour of the molecular 
correlation functions for 1,1' ^ 0. The correlator Sn(q) 
exhibits a pronounced maximum at q « 2.8, a q— value 
corresponding to about twice the average nearest neigh- 
bour distance. 

A further point worth mentioning is that all diagonal 
correlators Su(q) with I > 2 show a maximum at q = 
as was already found for 622(9) by Kammerer et al. p6[ . 
This might indicate a tendency for a local nematic order. 

In conclusion of this section we want to point out that a 
comparison of our findings with an analysis of the system 
in real space would be very valuable for a better under- 
standing of the molecular correlation functions. It would 
be especially interesting to connect the shifts in peak po- 
sitions, the relation between the peak amplitudes of dif- 
ferent correlators and the peak shapes to a microscopic 
characterization of the system in real space. 



B. Neutron Scattering Function 

As already mentioned above the results of section y are 
best tested for the static case. The molecular correlation 
functions necessary for an evaluation of the series rep- 
resentation eq.(|5j) have been presented in the previous 
subsection for I, I' < l co — 4. The result of this calcu- 
lations will be compared to the exact data for neutron 
scattering as determined according to eq. (0) . 

Since we are dealing with the special case of linear 
molecules and static correlations eqs. (p|)— (|7|) can be fur- 
ther simplified. As justified in the previous subsection 
the indices n, n' can be set to for the coherent terms 
if the z-axis of the body-fixed frame of reference is cho- 
sen in direction of the symmetry axis connecting both 
atoms of the molecule. The polar coordinates of the 
atoms A and B are given by (r^, i9a, Pa) — (d/2, 0, 0) and 
(rB,$B,<fB) = (rf/2,7r,0) respectively. This information 
can be inserted into eq.(0) for the coherent molecular 
form factors. Taking into account the trivial expression 

&ln I'n'il) = $11' firm' f° r the static self correlations and us- 
ing sum rules for the spherical harmonics and spherical 
Bessel functions p3 all summations over the incoherent 



terms can be carried out leaving only one q-dependent 
function b mc (q) usually termed molecular structure fac- 
tor. We can combine these reformulations to get the 
following result for the ns-static structure factor of di- 
atomic rigid molecules. 



S n (q)^St°Hq)= £ bW h (q)S$(q) + b™ c (q) 



U'<l. 



b% h (q) = [(2l + l)(2l' + l)}ij l (^)j l ,(^) ■ 
■ \ (<'h + (-1)'^) (S + (-1)''^) 



(9) 

(10) 

(11) 



+ a coh a cohMqd) 

We just want to mention that choosing l co = yields 
the usual Sears-expression pf. This fact demonstrates 
again that the representation through molecular correla- 
tion functions is a natural extension of the approach by 
Sears. 

The molecular structure factor b mc {q) shown in Fig. 
for the system of diatomic molecules with a special choice 
of scattering lengths has a very simple q-dependence. 
Starting from a maximum at q = it quickly decays 
and shows oscillations around the asymptotic value for 
q — > oo which is given by the g-independent first term 
of eq.(jll|). In the case of general molecules the q— 
dependence is given by a linear superposition of functions 
3q{qTw') involving all intramolecular distances r vv i = 

Both, the molecular correlation functions Sw (q) and 
the coherent form factors bf^, {q) are symmetric with re- 
spect to I and V . Thus the summation in eq.(|9|) can be 
restricted to V > I. If the scattering lengths of A- and 
B-atoms are equal the coherent molecular form factors 
b1i, h (q) vanish unless I and V are both even, i.e. only 
molecular correlation functions with I and V even con- 
tribute to S n (q). This is a specific property of the di- 
atomic linear molecule closely related to head-tail sym- 
metry. In case of exact head-tail symmetry the correla- 
tors S^, (q) are zero for I or I' odd. If only the scattering 
lengths are equal such correlations may exist but they do 
not enter S n (q) since the neutron probe can not distin- 
guish between A- and B-atoms. 

Apart from this system-specific property eq.dK) and 
the examples of form factors presented in Fig. [3 clar- 
ify some properties of the g-dependence of bf°, (q) which 
are quite universal, i.e. valid also in the time-dependent 
case and for general molecules. In both subfigures one 
can observe that the weight 6^°' l (q) of a given correlation 
function Su> (q) depends very strongly on the value of q. 
For instance the center of mass correlations 5*00(9) are 
very important for q »s but only of minor relevance for 
q > 10. From the mathematical properties of the spher- 
ical Bessel functions, determining the ^-dependence of 



b1°, h (q) also in the general case given by eq.(M), one can 
conclude a systematic relation between the values of I 
and V and the range of g-values in which the correla- 
tor Sw(q) will have a large weight. This relations are 
demonstrated in Fig. ||. Figure (a) shows the diagonal 
form factors 6^ o/l (g) up to I = 4. All of them are pos- 
itive over the entire g-range. bQQ h (q) has a maximum 
for q = which decays to zero with increasing q and 
is followed by further oscillations with strongly reduced 
amplitude. The other diagonal correlators also show a 
pronounced maximum that is followed by smaller oscilla- 
tions but the location of this main maximum is shifted to 
higher values of q. This shift grows monotonously with 
increasing I. Further we see that the maximum's ampli- 
tude decreases and its width gets larger upon increasing 



I. We have chosen a^ oh 



1.4 different from af oh = 0.25 



in order to have nonvanishing form factors for / or I' odd. 
Still we can observe that the dependence of bf°, (q) on the 
scattering lengths leads to smaller amplitudes for odd I. 
Fig. g (b) also shows the systematic behaviour for the 
offdiagonal form factors bQ°! l (q) with I = 0. Upon grow- 
ing V — I the amplitude of the first maximum quickly 
reduces and the contribution at the second maximum 
gets more important. Thus while the offdiagonal terms 
can not be neglected the relevant range is also shifted to 
higher values of q for growing V . The same is true for the 
nondiagonal terms with I > 0. 

Consequently the systematic behaviour of the weights 
bf°, h (q) can be formulated as a rule of thumb: The larger 
the values of I and I' — I the higher the value of q at which 
the form factor bf° t (q) will become relevant. We merely 
note that for the form factors as well as for the molec- 
ular correlation functions themselves negative values are 
possible for the offdiagonal terms. 

Having discussed all the relevant terms for the evalua- 
tion of the right hand side of eq. (B) we can now turn to 
the ns-static structure factor itself and test the conver- 
gence of the series representation, i.e. the quality of the 



different approximations Sn e °'(q), I. 



0,1, 



Fig. H (a) shows the exact result for S n (q) as eval- 
uated according to eq.(|l|). The scattering length have 
been chosen as af oh = 1.4, af oh — 0.25 and a" nc = 0, 
v = A, B. At q m 3 the static structure factor exhibits a 
small but well pronounced prepeak followed by a strong 
maximum at q « 6.6 and further maxima at q rs 13, 
q w 18 and q « 25. Thus the g-dependence of S n (q) 
shows the general structure also found for the molecular 
correlation functions. Taking a closer look to the various 
maxima and taking into account the properties we found 
for the weights bi°, h (q) we can further illuminate their 
origin. In case of the prepeak this is rather obvious since 
i) at q « 3 only correlators with small I and /' have to 
be taken into account (see discussion above) and ii) only 
Su{q) possesses a significant peak at q « 3. Therefore its 
origin must be connected to the structure found in the 



correlation function 5n(g). A quantitative analysis of 
this statement will be given in the following subsection. 
From the Figure || we can see that also at the position 
q rj 6.6 of the main peak correlators with I or /' larger 
than 2 do not play a big role. Thus the center of mass 
correlation function 5oo (g) and the molecular correlation 
function So2(q),S22(q) as well as 5oi(g) can be identified 
as the origin of the main peak in S n (q). The latter cor- 
relators being responsible for the slight shift of the peak 
position to higher q in comparison with the center of 
mass correlations. For the peaks at q rj 13 and q « 18 a 
number of correlators will be relevant. Using the weights 
tff, h (q) one can still conclude that S22(q) will be one of 
the main sources for the peak at q « 13 and S24 (g) as 
well as 544(g) will be of great relevance for the structure 
at q « 18. From the molecular correlation functions we 
have determined 544(g) is the one which contributes most 
to the peak at q sa 25 though correlation functions with 
1,1' > 4 could also be of great importance. 

Apart from the exact result for S n (q) Figure [| (a) also 
shows the result of the weak hindering approximation as 
used by Sears M . Except for the limiting values for q — > 
and q — > 00 the results of this approximation are rather 
poor since almost all of the structure in S n (q) is missed 
and even the main peak is not reproduced properly. This 
finding strongly supports results of Chen et al. |12] who 
observed that a factorization of translational and rota- 
tional correlations is not suitable for the description of a 
supercooled liquid. 

Figure (b) shows the convergence of the series rep- 
resentation by presenting the dependence of the error 
S n (q) — Sn°(q) on l co . The large error of the Sears 
result corresponding to l co = is already reduced to 
statistical fluctuations in a range of q up to q RJ 9 by 
choosing l co — 2. For l co = 4 the static structure factor 
S n (q) is almost perfectly reproduced for the entire range 
of q- values at least up to q « 20. Since the g-range cho- 
sen here compares well to the values accessible in a real 
neutron scattering experiment we can conclude that the 
convergence of the series representation (@) is fast enough 
to make it attractive for practical purposes. The cutoff 
of the infinite sum at a relatively small value of l co is well 
justified, at least for our choice of model. The actual 
value of l co that is sufficient will depend on the size and 
form of the molecule as well as their interactions |18| . In 
accordance with the given rule of thumb increasing the 
value of l co merely corresponds to an enlargement of the 
range in which 5 n (g) is reproduced. 



C. Prepeak 

We will now turn to a special feature of the static struc- 
ture factor. As already mentioned in the last subsection 
S n (q) shows a prepeak at g w 3. Given that prepeaks or 
first sharp diffraction peaks (FSDP) have been studied 



mostly in connection with strong glasses like e.g. S1O2 
p0[ or ionic glassformers |19| this observation is quite sur- 
prising. In the system of diatomic rigid molecules studied 
here we neither find a network structure nor do we have 
long range electrostatic interactions. Since the proposed 
explanations for the prepeak usually rely heavily on the 
network structure and/or the effects of the ionic charges 
they can not be applied to the prepeak in the present sys- 
tem. Then what could be the mechanism leading to in- 
termediate range order in a system without an extended 
network or long range interactions? 

In the preceding sections we have pointed out the con- 
nection between the g-dependence of 5„(g) and the con- 
tribution of various molecular correlation functions. To 
understand the mechanism responsible for the formation 
of a prepeak it will surely be helpful to analyze in detail 
which of the relevant molecular correlations contribute at 
q pp rj 3. The contributions of the intramolecular struc- 
ture factor b mc {q) and the center of mass correlations 
5q (g) are of large amplitude but they contain no struc- 
ture at all at q pp . Since they have opposite signs they 
almost cancel each other. Together they give the Sears- 
result which does not exhibit any prepeak (cf. Fig.[| (a)). 
The other relevant contributions at q pp are shown in Fig. 
H Thus the importance of the correlator 5n(g) for the 
formation of the prepeak is corroborated by this quan- 
titative analysis. Besides the negative contribution of 
5oi (g) no other term shows any structure at q pp . It may 
be noted that the negative contributions of 5n(g) and 
5o2 (g) at g rs 5 induce a better separation between the 
prepeak and the main peak. In accord with our "rule of 
thumb" the terms involving 1,1' > 2 are even less impor- 
tant than 522(g)- 

The predominant role of 5n(g) and the fact that no 
structure is observed in the center of mass correlation 
strongly suggest that orientational effects are responsi- 
ble for the prepeak in the studied system. Apart from 
the evidence given here it seems not too far fetched to 
consider spatial ordering due to sterical hindrance as a 
mechanism for prepeaks in molecular systems as close 
packing |3(]] and size effect |3l| have also been shown to 
be of great relevance in other systems. The importance 
of orientational degrees of freedom for the structure of 
5„(g) for liquid halogens and the prepeak in the molec- 
ular liquid CCI4 has also been put forward by Misawa 
p2Lp3f. Whereas in his work orientational correlations 
between neighbouring molecules had to be introduced as 
an assumption we have been able to give direct evidence 
for their relevance. 

Taking into account the interpretation of the molecular 
correlation functions (see above) it is also possible to get 
additional information on the kind of intermediate range 
order. 2?q cx Y® has the shape of a dumbbell which 
is, due to the choice of the body fixed frame, aligned to 
the symmetry axis connecting the atoms of the molecule. 



Thus the correlations of Sn(q) at q pp can be attributed 
to a preference for a parallel orientation of next nearest 
neighbours. Apart from the geometric effects leading to a 
parallel alignment of molecules also energetic affects are 
important for this " antiferromagnetic" order 



ABA 

I I I 
BAB 



(12) 



since the vicinity of A and B molecules is favoured by 
the choice of Lennard-Jones energies €ab < £aa — £bb ■ 

Prepeaks have often been reported to show special 
behaviour. In many glassy systems an increase of the 
prepeak amplitude with increasing temperature is found 
pOJplfl which is quite contrary to the behaviour of the 
other peaks in S n (q). The amplitude and position was 
also shown to be affected by the pressure [g0[ and the 
composition of the liquid |19[ . Therefore it will be inter- 
esting to characterize further the behaviour of the pre- 
peak in our system. 

Fig. |5| shows the temperature dependence of the pre- 
peak in the range T = 0.48 to T = 0.70. To eliminate 
the background given by the low q wing of the main peak 
we have substracted the S n (q) at T — 0.85. At this tem- 
perature the prepeak is not discernible anymore. With 
increasing temperature the prepeak decreases in ampli- 
tude and the width at half maximum grows. In the given 
temperature range no shift in the position of the peak can 
be observed whereas the main peak is shifted downwards 
by Aq = 0.2. Although seemingly uncommon for a pre- 
peak such temperature dependence has been reported be- 
fore. It has been found for carbon tetrachloride p4] and 
also in a recent study of m-toluol and m-fluoranilin by 
Morineau et al. 122] . Like our system these are molecular 
liquids. Thus one is tempted to infer that this behaviour 
might be typical for such systems. 

We are not able to determine the pressure dependence 
of the prepeak or to study effects of composition, i.e. the 
influence of atomic size and interactions since this has 
not been done in the simulation p6| . Still we want to 
point out one aspect which is closely related to compo- 
sitional studies but has never been considered before. A 
substitution of atoms of a certain species in the liquid 
will not only alter the interactions and the size of the 
atoms but will also affect the scattering lengths. Fig- 
ure M shows that the observability of the intermediate 
range order as a prepeak can indeed be influenced. As 
can also be inferred from eq.(||) in the case a^ oh = af oh , 
i.e. x = acoh/ a coh = 1) tne correlator Su(q) gives no 
contribution to S n (q) and the prepeak vanishes. Keep- 
ing the overall normalization constant and increasing the 
asymmetry in the scattering lengths the prepeak grows 
to its maximum value. Since at the small g-value of the 
prepeak position the molecular correlation function with 
I = 1 will be the most important one (apart from Soo(q)) 



a prepeak in molecular systems will be most easily ob- 
served for molecules with no head-tail-symmetry. Thus 
it is for instance not surprising that in liquid halogens no 
prepeak is observed although their liquid structure might 
be not too different from our system B2]. 



IV. CONCLUSION 

Tackling the problem of neutron scattering in super- 
cooled molecular liquids we have shown that well known 
concepts can be generalized naturely to give a concise 
description including all molecular degrees of freedom. 
Thus it was possible to reduce the large discrepancy that 
was found between the exact result for S n (q) and a de- 
scription in terms of the center of mass only. The result- 
ing description of S n (q,t) in terms of molecular corre- 
lation functions SJ^ i' n '(q,t) was shown to shed light on 
the (/-dependence offering clear relations between the q- 
range and the angular resolution described by the " quan- 
tum numbers" / and n. Giving a relation between the 
intermediate scattering function and the molecular cor- 
relations the representation could be used to extract (par- 
tial) information on the correlations of TDOF as well as 
ODOF and their interference from neutron scattering ex- 
periments. From the general expression we were able to 
give hints for an effort in this direction. A quantitative 
test of the formalism for a liquid of diatomic molecules 
in the "worst case" of static correlations led to a good 
agreement with exact results for S n (q) if molecular cor- 
relation functions up to l co = 4 are taken into account. 
This analysis also confirmed that a factorization like the 
weak hindering approximation which leads to a descrip- 
tion in terms of the center of mass only is not suitable in 
the supercooled regime. 

In the simulation of a liquid of diatomic Lennard-Jones 
molecules we observe a prepeak in the ns-static structure 
factor S n (q) which could be attributed to intermediate 
range orientational order. The temperature dependence 
of this prepeak is in accord with the results found for the 
molecular liquids CCI4, m-toluidine and m-fluoroanilinc 
but in variance with the behaviour found for most 00- 
valent glassformers. The influence of scattering lengths 
on the observability of intermediate range (orientational) 
order was examined offering the conclusion that a mani- 
festation as a prepeak may not occur in case of head-tail 
symmetry of the molecules. 
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FIG. 1. Molecular static structure factors Sut(q) for a system of diatomic rigid molecules at the lowest temperature T=0.477 
of the simulation. 
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